Graphene in a Strong Magnetic Field: Massless Dirac Particles vs. Skyrmions 
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We discuss models for massless Dirac fermions being subject to a perpendicular magnetic field 
in spherical geometry. These models are analogues of Haldane's spherical construction for massful 
charge carriers. The single particle states constructed here are easily implemented in existing numer- 
ical code for many-body problems in conventional quantum Hall systems. Moreover, the many-body 
states of fully filled sublevels in the subspace of lowest Landau level index are skyrmions with respect 
to the layer spin. 
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I. INTRODUCTION 

In the recent years, graphene has developed to one of 
clearly most active directions of work in today's both ex- 
perimental and theoretical condensed matter physicsi'^'^. 
Compared to conventional two-dimensional electronic 
systems, the peculiar properties of graphene mainly stem 
from its linear dispersion near the Fermi energy, and 
the chiral nature of electronic states entangling the mo- 
mentum and sublattice degree of freedomSi'^. Among a 
plethora of interesting phenomena, the quantum Hall ef- 
fect occurring at anomalous ("half-integer") filling fac- 
tors is one of the most spectacular observations in this 
new type of material^i^i^. Other partially related unusual 
features of graphene include a cyclotron mass being pro- 
portional to the square root of the particle density^, and 
a non-equidistant Landau level spectrum proportional to 
the square root of the magnetic field^i^i^ii^i^. 

In quantum Hall physics, the exact numerical treat- 
ment of finite many-body systems has always been an 
important source of theoretical evidenceiS. A particu- 
larly convenient model for such numerical simulations 
was given already 25 years ago by Haldan&ii. In this 
construction, massful electrons move on a spherical sur- 
face in a radial monopole magnetic field. This model has 
the advantage of lacking any system edge (and therefore 
reducing finite-size effects in numerical results). On the 
other hand, the mathematical properties of the states 
in the lowest Landau level are particularly simple, en- 
abling also substantial analytical progress. Both as- 
pects have made the Haldane sphere to a widely used 
model in theoretical, in particular numerical descrip- 
tions of quantum Hall systems; for some representative 
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Moreover, most recently this spherical model, explicitly 
constructed for massful electrons in conventional two- 
dimensional systems, was also applied to the situation 
of massless carriers in grapheno^2i^i2^i^i2i. In particu- 
lar, the authors of Refj^ state that "the analogous solu- 
tion for carriers with linear dispersion ... is not known 
for the spherical geometry" . The purpose of this note 
is introduce such models of massless charge carriers in 
the spherical geometry. As we shall see below, they are 
rather straightforwardly developed from the situation of 



massful particles. A particular feature of these mod- 
els are many-body ground states having the features of 
skyrmions known from conventional quantum Hall mono- 
layers. Finally we mention previous work by other au- 
thors on Dirac fermions coupled to gauge fields in spehri- 
cal geometry inspired by fuUerrenes^^'^'^, i.e. another 
allotrope of carbon . Differently from the models to be 
discussed here, these constructions do not rely on angu- 
lar momentum operators and have therefore also different 
spectra. 

This paper is organized as follows. In section [TT] we re- 
view the basic properties of the spherical model for mas- 
sive charge carriers in a perpendicular magnetic field. In 
section IIIII we briefly recall elementary facts concerning 
planar graphene in a perpendicular magnetic field, be- 
fore constructing analogous models in spherical geometry. 
Here we also discuss two-body and many-body states. 
We close with conclusions in section HVl 



II. 



SPHERICAL MODEL FOR MASSFUL 
CHARGE CARRIERS 



In Haldane's spherical model for massful charge carri- 
ers in a perpendicular magnetic field, electrons move on a 
spherical surface of radius R which is penetrated by a ra- 
dial monopole magnetic field B = hcS/eB? . Here e > 
is the elementary charge, and 2S is the integer number 
of flux quanta hc/e through the surfacoii, i.e. the radius 
is given by i? = i^fS where i = \fhcfeB is the magnetic 
length. The single-particle Hamiltonian reads 
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Here M is an effective mass, uj = eB/Mc is the cyclotron 
frequency, and the kinetic angular momentum reads, us- 
ing again standard notation. 



A = r X 



P 



-A 



(2) 



where the vector potential yields V x A = Bft, VL ^ r/R. 
The operator A has the elementary properties 



A-fi = Sl-A = 0, 



(3) 
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[A", A^] = ihe"'^'^ (A'^ - hsn"') . 



(4) 



In particular, the components of A fail to fulfill proper an- 
gular momentum commutation relations and can there- 
fore not be considered as the generators of rotations. In- 
stead, rotations are generated by the angular momentum 
operator 



fulfilling 



L = A + hsn 



L-n = n- L = hs, 



= + {hsy 



(5) 



(6) 
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Thus, the eigenvalues of are given by h^l{l + 1) with 
I = S + n and n S {0, 1, 2, ...} where n = corresponds 
to the lowest Landau level. Note, however, that the spec- 
trum of the Hamiltonian ([1]) is not equidistant, but grows 
quadratically with n, differently from the planar situa- 
tion. 

Let us choose the usual gauge A = {hS /eR)ip cot 'd, 
where are the usual polar coordinates and if) is the 
unit vector in the azimuthal directio n^^i^^ . Then eigen- 
states in the lowest Landau level have a particularly sim- 
ple structure given by 
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with TT = p + eA/c and v w lO^m/s. The Pauli matri- 
ces describe the sublattice or pseudospin degree of free- 
dom, and the Zeeman coupling to the physical electron 
spin has been neglected. The double sign (±) (valley 
index) determines which corner of the Brillouin zone in 
considered^^. These two cases are, in the absence of a 
magnetic field, related by time reversal^ implemented by 
a complex conjugation, i.e. changes sign while cr^ 
remains unaltered. Note that this behavior is different 
from angular momentum operator describing a proper 
spin and not a sublattice degree o freedom. In what fol- 
lows we shall concentrate on the case {+). Defining the 
usual bosonic operators 
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fulfilling [a, a^] = 1, the Hamiltonian reads 



n 



pi 
(+) 



hv r- 



aa 



a a 



(13) 



where = {a^ ± i(T^)/2. The well-known eigenstates^ 
of the Hamiltonian are given by |0,t) with energy 
£0 = and, for n > 0, 



\n,±) 



1 
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(|n,T)±|n-l,i)) 



(14) 



with energy = ±{hv/£)\/2n. Here n is again the Lan- 
dau level index, and the arrows are obvious standard no- 
tation for the sublattice spin states. In particular. Lan- 
dau level index and (sublattice) spin are entangled in 
these eigenstates, a feature that will be reproduced by 
the spherical models for graphene in a magnetic field to 
be discussed now. 



where u = cos(^9/2)e''^/^ v = sin(i?/2)e-*'^/2^ ^nd m e 
{— .S, . . . 5} is the eigenvalue of /h. 

We note that the gauge invariant angular momentum 
operator ([5]) crucially depends on the magnetic fiux and 
is not identically to the gauge-dependent canonical angu- 
lar momentum operator Lean = r x p. Moreover, due to 
the relation ([5]), one could alternatively define the Hamil- 
tonian 



, 1 L2 

which differs from ([T]) just by a trivial constant. 



(10) 



A. Spherical models 

Let us now discuss models for massless Dirac particles 
on a sphere penetrated by a radial monopole magnetic 
field, for which we shall use the same gauge as stated 
before. Note that the crucial step in constructing the 
spherical model for massful carriers is to replace the lin- 
ear momentum of the usual planar Landau problem with 
an appropriately defined angular momentum. This shall 
also be our guideline regarding the case of massless Dirac 
fermions. 



III. MASSLESS DIRAC PARTICLES IN 
SPHERICAL GEOMETRY 

For a planar graphene sheet in a perpendicular mag- 
netic field, the single-particle states around one of one 
of the two inequivalent corners of the first Brillouin zone 
are described by^ 



(±) 



(11) 



1. Single-particle Hamiltonians 

Inspired by the expressions ([1]) and (fTO)) we consider 
the Hamiltonians 



and 



^(±) = . f(±) = ^A . r(±) 



(15) 
(16) 
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Here the operators describing the sublattice degree 
of freedom are again just given by the PauH matrices, 
f(+) = (7, while r(_) are the negatives of their complex 
conjugates, T(_) = —((?)*, in close analogy to the planar 
case ([TT|). Note that both sets of operators, although not 
referring to a proper spin, fulfill the usual relations 



(±) 



(17) 



This fact will chiefly facilitate the analysis of the Hamil- 
tonian (|16|) in terms of elementary angular momentum 
theory. The treatment of the Hamiltonian p^ . however, 
is more complicated since the components of A do not 
fulfill an angular momentum algebra. Moreover, in the 
following we shall again concentrate on the case (+), the 
case (— ) can be treated completely analogously and is 
just related via complex conjugation. 

The two Hamiltonians given in Eqs. p^ . ([15]) differ 
by the radial contribution 



(+) 



hv 



(+) 



(18) 



spin 1/2, one can explicitly formulate the 25 + 2 eigen- 
states in the multiplet j = I + 1/2, 



S+l/2 + m 
25 + 1 

5+ 1/2-m 
25 + 1 



'Pm-i/2{u,v)\ t) 

f,n+l/2{u,v)\ i) , (21) 



where m e {—S — 1/2, . . . , 5 + 1/2} is the eigenvalue of 
J^' /h, and the 25 states with j = I — 1/2 read 



V'm(",w) 



5+1/2-J71 
25+1 



<Pm-l/2{u,v)\ t) 



5+ 1/2 + m 
25 + 1 



Vrn+1/2{U,V)\1) (22) 



with me {-5+1/2, 5-1/2}. Similarly to the planar 
model of graphene, the sublattice spin and the conven- 
tional orbital motion are entangled with each other. 

Denoting ^^{u,v) = (r|m,±), the expectation values 
of i7 • (7 are straightforwardly calculated as 



which is, differently from the expressions (fTO|) and ([T]), 
not just a trivial constant. In fact, one might argue 
that should be considered to be closer to the planar 

model of graphene since the operator A, in contrast with 
L, does not have a radial component. However, as we 
shall see shortly, when concentrating on the subspace of 
lowest Landau level index n = 0, the solutions of nega- 
tive energy are simultaneous eigenstates oi L- a and il ■ a 
(and therefore also A • (?) at any systems size, whereas 
the solutions of positive energy become such simulta- 
neous eigenstates in the thermodynamic limit S oo. 
Thus, as far as the subspace of lowest Landau level in- 
dex is concerned, the eigenstates of and '^(+) arc 
either identical at any system size or become identical 
in the thermodynamic limit, where the planar model of 
graphene is recovered. Let us therefore first concentrate 
on the latter Hamiltonian. 

Introducing a total angular momentum of the usual 
form J = L + ha/2 one can rewrite the Hamiltonian m 
as 



(+) 



(19) 



Moreover, J commutes with the Hamiltonian and admits 
total angular momentum quantum numbers j = l±l/2 — 
5 + n ± 1/2. Thus, the spectrum of Ti.',,-. reads 



± — p= 5 + ?i. + - + - 

eVs\ 2^2 



(20) 



Let us again focus on the lowest Landau level index n = 
0. Using the well-known Clebsch-Gordon coefficients^^ 
for coupling an angular momentum of length 5 with a 



(m, ■ (T|m, +) 
(m, —\Q ■ aim, — ) 



25 
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-1 . 



(23) 
(24) 



Thus, the variance of Q ■ a within the states of lower 
energy \m, —) is exactly zero, while for the states |m, +) 
one finds 



rn, - 



1 - 



(25)^ 



(m, +\n ■ (T\m, +) 



25+1 



(25 + 2)^ (25 + 2)^ 



(25) 



As a result, the eigenstates |m, — ) of the Hamiltonian 
7i'^^^ are for any system size simultaneously also eigen- 
states of whereas the eigenstates |m,+) achieve 
this property in the thermodynamic limit 5 — > oo. In 
particular the energy e_ of corresponding to e'^ 

for n = is just e_ = —hv/i\/S. Moreover, replacing, as 
an approximation, the operator f2 • a with ±1 one finds 
for the energy spectrum of 



hv 



1 1 

2^2 



(26) 



This result should be seen in analogy to the conventional 
spherical model ([T]). Here the energies increase quadrat- 
ically with the Landau level index n, while the spectrum 
of the underlying planar model is of course equidistant 
with energies proportional to n. In the case of graphene, 
the energies of the planar model are proportional to 
while in the spherical model 'H{+) they increase with 

^Jri' = n. As seen above, the approximation underly- 
ing Eq. is exact for both branches of the spectrum 
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at n = and S —f oo. It is an interesting specula- 
tion whether it becomes also exact in the thermodynamic 
limit if n ^ 0. 

The single-particle states given above were also exam- 
ined some time ago by Rezayi in circumstances of quan- 
tum Hall skyrmions^, an analogy we shall explore in 
some detail in section UlI A 31 



2. Two-body states and interaction matrix elements 

For the conventional spherical model for massful charge 
carriers described in section|TTl a particularly simple form 
for two-body states of given total angular momentum can 
be devisedii— . Using these expressions, matrix elements 
of rotationally invariant interactions are conveniently pa- 
rameterized in terms of pseudopotentialsi^. 

In the present case of massless particles, the pseu- 
dospin degree of freedom adds to the complexity of the 
two-body problem, and we have not found a similarly 
concise expression for states with good quantum numbers 
of the total angular momentum. However, what is usu- 
ally needed in numerical implementations of interaction 
Hamiltonians are matrix elements between tensor prod- 
ucts of single-particle states which are related to pseu- 
dopotentials via Clebsch-Gordan coefficientsi^. Now, us- 
ing the expansions pT|) . ((22l) it is straight forward to ex- 
press such interaction matrix elements of states of mass- 
less particles considered here in terms of those of massful 
objects given by Eq. Thus, regarding numerical im- 
plementations of interaction operators, it is an easy and 
straightforward task to adjust an existing code for the 
conventional spherical model to the single-particle states 
of massless carriers. 



3. Many-body states 

We now discuss many-body states of fully filled ener- 
getic sublevels (PU|) of lowest Landau level index = 0. 
Let 1^^) denote the Slater determinant of all 25 -I- 1 ± 1 
single-particle states given in Eq. pT|) and ([22]) . respec- 
tively. Both many-body states are singlets of the total 
angular momentum. For the particle density one natu- 
rally finds 



n±(f) = (*±|^5(f-f,)|*±) 

i 

2S'+l±l 
while the pseudospin density is given by 



(27) 



2^ / sin'd cos If 



(28) 



in accordance with Eqs. ([23)) . ((24|) . We see that these 
pseudospin densities form the typical "hedgehog" struc- 
tures known from skyrmions^Siii. However, differently 
from skyrmions in conventional quantum Hall monolay- 
ers, the physical electron spin as well as the valley spin 
are polarized here while the sublattice spin is forming 
a topologically nontrivial structure. In fact, the single 
particle wave functions (PT|) . along with the many- 
body states where discussed already in Ref3 as 
models for (electron spin) skyrmions. Here we have iden- 
tified formally the same skyrmion states (with respect to 
the sublattice spin) as ground states of fully filled Lan- 
dau levels of massless Dirac particles. We stress the fact 
that the Hilbert spaces spanned by the single-particle 
states ([2l]) . ([22]) for massless Dirac particles are differ- 
ent form the Hilbert space spanned by the single-particle 
wave functions ([9]) for conventional massful carriers. In 
several recent publications however, the Landau levels of 
graphene were modeled by the conventional wave func- 
tions for massive particles^Sii^'^i^^'^. 

Finally, in order to evaluate interaction terms within 
the many-body states |\E'''=) it is useful to consider the 
pair distribution function 

9^i\ri - r2\) = E '^(^^1 - ^^^)'^(^^2 - r;)|«'±) . (29) 



Here we obtain 

.9+(0 = 



+ , , _ (25 + 2)^ 



(477^25.)^ 



2 \ 2S-1 



1 - 



(25 + 2)2 V 2^225 
45 + 2 r2 



and 



9 (r) 



(25)2 



(471^2^)- 



(25 + 2) 



1-1- 



2£225 



r2 \ 
2£225 J 



2S-1' 



(30) 



(31) 



From these expressions the ground state energies for ar- 
bitrary two-body interactions can be evaluated by inte- 
gration. In particular, for Coulomb interaction we find 



E+ = - 



1 



14S-1 



V5 



45 
25 

24S-1 



5- 



45- 



45 + 1 



45 
25 



5. 



(32) 



(33) 



Here e is the dielectric constant of the host material, and 
we have as usual assumed that the direct (Hartree-) con- 
tribution to the grond state energy is cancelled against a 
neutralizing background. The result p3|) agrees with the 
one given ini^, while Eq. ([5^ differs in detail form state- 
ments made there by a contribution which, however, van- 
ishes in the thermodynamic limit. If the states l^*^) are 
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interpreted as charged excitations of a conventional quan- 
tum Hall ground state at a filling factor of unity, their 
energies should be compared with conventional ground 
state energy given by2& 



e^2S + l 



^4S 



45 
25 



5. 



(34) 



To provide a meaning full comparisoi*ii, one also has to 
take into account that the particle number in the excited 
states 1 4'*) differs from the conventional ground state by 
±1. Then one finds for the excitation gap in the thermo- 
dynamic limit 5 — !■ oo 

, 25+l±l e^l 
^ 25+1 ^ ^7Z4V2' 

in accordance again with Ref4S, which is exactly the 
result predicted by field-theoretical considerationsii. 
Moreover, in Ref.— it was also found numerically, that 
the state l^""*") has a vanishing variance for Coulomb in- 
teraction operator in the thermodynamic limit, i.e. this 
state becomes an eigenstate of the interaction operator in 
the limit of an infinite system. This surprising result was 
also reported for other types of long-range interactions^. 



Hamiltonian (fT5l) is arguably closer to the planer model 
of graphene, the latter one ([TB)) can be analyzed easily 
by elementary angular momentum theory. Both Hamil- 
tonians differ by a nontrivial operator. However, in the 
subspacc of lowest Landau level index n = the eigen- 
states of the single-particle Hamiltonian (|16p become also 
cigenstatcs of ([TS1) . The latter result holds for the states 
of positive energy in the thermodynamic limit, while for 
the states of negative energy this statement is true at 
arbitrary system size. In particular, the Hilbert spaces 
spanned by the single-particle eigenstates of the Hamil- 
tonians introduced here are different form the Hilbert 
space spanned by the single-particle wave functions for 
conventional massful carriers. It is a very interesting 
question for future work, whether recently reported re- 
sults of numerical studies of quantum Hall physics in 
graphene using conventional wave functions for massive 
particles^i^'^i^^'^ are possibly altered if single-particle 
states for massless charge carriers are used. Indeed, the 
single particle states constructed here are easily imple- 
mented in existing numerical code for many-body prob- 
lems in conventional quantum Hall systems. Finally, the 
many-body states of fully filled sublevels in the subspace 
of lowest Landau level index are skyrmions with respect 
to the sublattice spin. 
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IV. CONCLUSIONS AND OUTLOOK 

We have introduced spherical models for the mass- 
less Dirac charge carriers of graphene being subject to 
a perpendicular magnetic field. The Hamiltonians psp . 
(dni presented here are analogues of Haldane's spherical 
construction for massful charge carriers. While the first 
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